Quantum simulation of time-dependent Hamiltonians and the convenient illusion of 

Hilbert space 
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We consider the manifold of all quantum many-body states that can be generated by arbitrary 
time-dependent local Hamiltonians in a time that scales polynomially in the system size, and show 
that it occupies an exponentially small volume in Hilbert space. This implies that the overwhelming 
majority of states in Hilbert space are not physical as they can only be produced after an expo- 
nentially long time. We establish this fact by making use of a time-dependent generalization of the 
Suzuki- Trotter expansion, followed by a counting argument. This also demonstrates that a com- 
putational model based on arbitrarily rapidly changing Hamiltonians is no more powerful than the 
standard quantum circuit model. 



The Hilbert space of a quantum system is big — its di- 
mension grows exponentially with the number of parti- 
cles it contains. Thus, parametrizing a generic quantum 
state of TV particles requires an exponential number of 
real parameters. Fortunately, the states of many physical 
systems of interest appear to occupy a tiny sub-manifold 
of this gigantic space. Indeed, the essential physical fea- 
tures of many systems can be explained by variational 
states specified with a small number of parameters. Well 
known examples include the BCS state for superconduc- 
tivity [l|, Laughlin's state for fractional quantum Hall 
liquids 2\, tensor network states occurring in real-space 
renormalization methods j3j|. In these cases, the number 
of parameters scales only polynomially with N. 

In this Letter, we attempt to define the class of phys- 
ical states of a many-body quantum system with local 
Hilbert spaces of bounded dimensions, and prove that 
they represent an exponentially small sub-manifold of the 
Hilbert space. We say that a state is physical if it can be 
reached, starting in some fiducial state (e.g. a ferromag- 
netic state, or the vacuum), by an evolution generated 
by any time-dependent quantum many-body Hammilto- 
nian, with the constraint that 1) the Hamiltonian is local 
in the sense that it is the sum of terms each acting on 
at most k bodies for some constant k independent of N, 
and 2) the duration of the evolution scales at most as a 
polynomial in the number of particles in the system. The 
assumption about the initial fiducial state is artificial; we 
could alternatively define the class of physical evolutions 
for quantum many-body systems as the ones generated 
by Hamiltonians obeying constraints 1 and 2, and would 
reach the same conclusions. 

The second constraint is very much reminiscent of the 
way complexity classes are defined in theoretical com- 
puter science, where that central object of study is the 
scaling of the time required to solve a problem as a func- 
tion of its input size. The classical analogue for the prob- 
lem that we address is a well known counting argument 
of Shannon ^ demonstrating that the number of boolean 



functions of N bits scales doubly exponentially (as 2^ ), 
with the consequence that no efficient (i.e. polynomial) 
algorithm can exist to compute the overwhelming major- 
ity of those functions. Indeed, the number of different 
functions that can be encoded by all classical circuits of 
polynomial depth scale as 2^°'*''^^' , which is exponentially 
smaller than the total number of Boolean functions. 

Our contribution is a quantum generalization of this 
result, which has in part been a folklore theorem in the 
quantum information community for some time. The 
crux of our argument is to demonstrate that the dy- 
namics generated by any local Hamiltonian, without any 
assumptions on its time-dependence, can be simulated 
by a quantum circuit of polynomial size. All previously 
known simulation methods (5l-[To| produced a quantum 
circuit of complexity that depends on the smoothness of 
the Haniioltonian, scaling e.g. with \\dH/dt\\ or some 
higher derivatives. Using the results of Huyghebaert and 
De Raedt ll| , we show how this conditions can be over- 
come, and demonstrate how the counting argument of 
Shannon can be repeated in the quantum case by invok- 
ing the Solovay-Kitaev theorem [12;]. Note that a direct 
parameter counting would not produce this result be- 
cause we impose no restriction on the time-dependence of 
the Hamiltonian. The complete description of a rapidly 
changing Hamiltonian requires lots of information, so 
from this perspective there are in principle enough pa- 
rameters to reach all states in the Hilbert space. This 
leads to the conclusion that most states in the Hilbert 
space are not physical: they can only be reached after 
an exponentially long time. This has to be contrasted to 
the classical case, where all states of N bits correspond 
to physical states: they can easily be generated by trivial 
depth-one circuits. The difference between classical and 
quantum behavior is due to the existence of quantum 
entanglement. 

Our demonstration that arbitrary local time- 
dependent Hamiltonians can be efficiently simulated on 
quantum computers is of interest in its own right in the 
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context of quantum computation. More precisely, we 
are concerned with Hamiltonians acting on N particles 
of the form 



Hit) 



(1) 



XC{1,2,...,N} 



formula to time-dependent Hamiltonians without com- 
promising the error, where the Hamiltonian may exhibit 
fluctuations much faster than the time step At. We begin 
by breaking the total time evolution into short segments 
t/(0, T)) = C/(t„, i„ + At)... U{t2M + At)U{Q, + At) 
duration of time At 



where X labels subsets of the N particles, each term 
has bounded norm ||7fx(i)|| < E, and each term acts on 
no more than k particles, i.e., Hx{t) = if \X\ > k, 
and k is fixed independent of the system size. We 
make no assumption on the geometry of the system and 
the coupling can be arbitrarily long ranged. The time- 
evolution operator U{0,t) from time to t is governed 
by Schrodinger's equation -^U{0,t) = ~iH{t)U{0,t), 
with solution given in terms of a time-ordered integral 

f/(0,t) =rexp{/oF(s)ds}. 

Starting with Feynman's exploration of quantum com- 
puters [Ml, it has been well established that the time- 
evolution operator generated by Hamiltonians of the form 
Eq. ([1]) can be decomposed into short quantum circuits, 
provided that H{t) varies slowly enough (Bl-fioj. In all 
cases, this is achieved by approximating the evolution 
operator by a product formula 



U{t 
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At) =T exp 



tj+At 



dsJ^Hxis) 

X / 



In the simple case where the sum over X contains only 
two terms, say Hi and H2, it has been shown that 
the generalized Trotter-Suzuki expansion 

C/^^ {tj , tj + A<) = r exp ( j dsHi (s) j 

/ rtj+At \ 

xTexp —i i dsH2{s) 

gives rise to an error in terms of operator norm that is 
equal to 



\\U{tj,tj+At) ~ U''='{tj,tj +St)\\ < ci2{Aty 



U{Q, i) « n e^v{-^Hx^{tp)Atp] , (2) ^ith ^12 of the order of 1 and given by 



where the sequences Xp, tp, and Atp are set by spe- 
cific approximation schemes, such as the Trotter formula 
[Til l or the Lie-Suzuki- Trotter formula [IB]- Because each 
term Hx acts on at most k particles, this last expres- 
sion represents a sequence of fc-particles gates. A stan- 
dard quantum circuit is obtained by decomposing each 
of these fc-body operator as a sequence of one- and two- 
qubit gates using the result of Solovay-Kitaev [H, ■ 

Perhaps the simplest example of a product formula de- 
composition of U{0,t) is given by 



C12 



tj+At 



dv 



du\\[Hi{u),H2iv) 



U{0,t) ^l[l[cxp{-iHx{jAt)At} 



(3) 



where the product over X can be carried in any given 
order. This decomposition makes use of two approxima- 
tions. First, the time-dependence of the H{t) is ignored 
on time-scales lower than At: the Hamiltonian is approx- 
imated by a piece-wise constant function taking the val- 
ues H{jAt) on the time-interval [{j—l)At, jAt]. Second, 
each matrix exponential is decomposed using the Trot- 
ter formula exp{-iH{t)At} « eyip{-iHx{t)AtHx}. 
Clearly, the size At of the time intervals must be shorter 
than the fluctuation time-scale of H{t) for the flrst ap- 
proximation to be valid. At ^ WdH/dtW"-^. Higher fre- 
quency fluctuations would therefore require breaking the 
time-evolution into shorter intervals, thus increasing the 
overall complexity of the simulation. 

Time- dependent Trotter- Suzuki expansion — Somewhat 
surprisingly, it is possible to generalize the Trotter-Suzuki 



which is upper bounded by with Cmax = 

maxt maxx Note that this bound does not de- 

pend on the derivative of the Hamiltonian (and is there- 
fore also valid for non-analytic time-dependence). Note 
also that the bound reduces to the usual Trotter error 
for the time-independent case and is therefore equally 
strong, and that it can straightforwardly be generalized 
to higher order decompositions. 

For our present application, the Hamiltonian is the 
sum of i e poly(A^) fc-body terms, c.f. Eq. ([1]). We 
can therefore iterate the above procedure log2(i) times; 
at the n'th iteration, there are 2" terms, each of strength 
upper bounded by CmaxL/2". The total error for ap- 
proximating the exact time-evolution of the Hamiltonian 
with L terms by a product of L time-ordered terms is 
therefore 



1 



.(At) 



m — 1 



L 



< icL.^'(At)^ 



which can be made arbitrary small by choosing a AT that 
scales as an inverse polynomial in the number of qubits. 
Approximating the time evolution operator over a total 
time T with a product of fc-local unitaries such that the 
total error is e/2 can therefore be achieved by choosing 



At: 



r 2 ■ 
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The total number G oi k -local quantum gates to achieve 
this accuracy is then equal to 



2 r3 



-T'L 



T c2 
^ ' ^ AT e 



Each of these /c-body unitaries can be decomposed into 
standard 2-qubit gates chosen from a discrete set of quan- 
tum gates (e.g. CNOT's between any pair of qubits sup- 
plemented by a local tt/S rotation gate) . The obvious 
approach is to exactly integrate the corresponding time- 
ordered exponential on a classical computer, and then 
to use the Solovay-Kitaev theorem to count the num- 
ber of standard quantum gates needed to approximate 
this. To achieve an accuracy ec per unitary, we need 
dsK (}og''^'^ {1/ ea)) standard gates with csk and dsK 
constants; we choose cq such that 



2G(e,T)- 



The total number of quantum gates as chosen from a 
discrete set of gates needed to approximate the complete 
time evolution with an error e is therefore upper bounded 
by 

Gtot{t,T) = dsKG{e,T)\og^''< (G(e,T)/e) 

which is polynomial in the number of qubits; it roughy 
scales quadratically in T and as the cube of the total 
number of non-commuting local terms in the Hamilto- 
nian. This actually proves that quantum computers op- 
erating under arbitrary time-dependent local Hamiltoni- 
ans are no more powerful than quantum computers based 
on the quantum circuit model. 

Average Hamiltonians and randomized evolution — Note 
that the time-dependent Trotter-Suzuki decomposition 
described in the previous section does not lead to 
a product formula because each term appearing in 
it involves a time-ordered integral Ux(tj,tj + At) = 

T exp{—i J*'^"^*^ Hx{s)ds} rather than the exponen- 
tial of a term of the Hamiltonian at a given time 
exp{— iAti?x(i)} as in Eq. Although this does 

not affect the conclusions reached in the next section 
on the counting of possible quantum states, it is un- 
satisfactory from the point of quantum simulation. In 
this section, we demonstrate how to recover a product 
formulae by making use of randomness. In Ref. [l7| . 
product formulae decompositions Eq. ([2]) were found 
for any local Hamiltonian, where the number of terms 
Np in the product depends on a smoothness parame- 

ter Ap = snpo<p<P,o<t<TExi\\dtHxmy^^''+''>- In 
particular, these decompositions are inefficient when the 
fluctuation time-scale of the Hamiltonian becomes too 
small. Our methods circumvents these requirements by 
using randomness. 

This is done in two steps. First, we can replace the 
time-ordered exponential integral with the exponential 



of an ordinary integral without introducing a significant 
error. Indeed, we show in Appendix \K\ that 



Texp < —i 



exp < — i / dsHx{s) 



<2\\HxfAt (4) 



Using this result, we obtain the approximate decomposi- 
tion 



C/(0,i) 



nn-p 

3 = 1 X 



tj+At 



Hx{s)ds 



(5) 



Note that this is still not a product formula because it 
involves integrals. This first step has nevertheless elimi- 
nated the need of a time-order operator. 

The second step to obtain a product formula for 
U{0,t) — one that does not require any integrals — makes 
use of randomness. The average Hamiltonian H^'^j on the 
interval [(j — l)At,jAt] can be estimated using Monte 
Carlo integration. For every j, we can pick m ran- 
dom times Tj e [tj , tj + At] and approximate « 
m SfeLi Hx{Tj)- Because the variance of the Hamilto- 
nian is bounded by H-ffjflP, the sum converges to -ff^j 
with error estimate At\\Hx\\/ \frn. Using this Monte 
Carlo average, we can approximate the evolution oper- 
ator of the time interval [tj , tj + At] by 

U^-(t,,t,+At)^exp{^t-J2Hx{T^)} (6) 

fc=i 

™ At 

k=l 

where the order of the product can be chosen according 
to increasing values of Tj. The error in the first ap- 
proximation Eq. (|6]) is set by the Monte Carlo estimate 
At\\Hx\\/y/m while the second approximation Eq. ([T]) is 
the usual Trotter-Suzuki. Summarizing, wc can decom- 
pose the total evolution operator from time to T as 



f/(0,T)« n eM-^—Hx{r^)} 

J,k,X 



(8) 



where the product should be taken in increasing order 
of Tj and any order of X. This is a standard product 
formula Eq. ^ — identical to the usual decomposition 
explained in the introduction Eq. ([3]) and the one pre- 
sented in — except that the times rj^ at which the 
Hamiltonian is sampled are random. Thus, we see that 
by sampling the Hamiltonian at random times, we com- 
pletely circumvent any smoothness requirements. 

We note that this proof of Eq. (jH) is an illustration 
of the decoupling principle that tells us that the high- 
frequency fluctuations of the Hamiltonian should not af- 
fect the low-energy physics. As a consequence, it is possi- 
ble to largely ignore these fluctuations — by replacing the 
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time-dependent Hamiltonian by its average value on each 
time bin — without significantly modifying the dynamics 
of the system. This is the working principle behind renor- 
malization group methods of quantum field theory and 
quantum many-body physics. The rotating wave approx- 
imation [l^ and effective Hamiltonian theory are 
simple examples illustrating this principle in the case of 
time-dependent Hamiltonians. 

More generally, we show in Appendix [B] that we 
can replace the time-dependent Hamiltonian H{t) with 
a smoothed version H{t) with fiuctuation time-scale 
bounded by a without significantly affecting the resulting 
time-evolution operator. More precisely, we show that 
the time-evolution operators from time to T differ by 
at most \\H\\^Ta. 

Counting states — Let us now consider the set of all quan- 
tum states that can be reached starting from some fidu- 
cial state |0) and evolving for some polynomial amount 
of time under any time-dependent Hamiltonian. A direct 
counting argument appears difficult because there are in- 
finitely many distinct time-dependent Hamiltonians, and 
therefore there can a priori be infinitely many states in 
that set. However, we have just established that the time 
evolution operator generated by any one of these Hamil- 
tonian can be well approximated by a polynomial-size 
quantum circuit. 

Quantum circuits are discrete objects. Because the 
Hamiltonian is fc-local, each term exp{—iAtHx{tj)} ap- 
pearing in the product formula Eq. is a fc-body quan- 
tum gate. The Solovay-Kitaev theorem shows that each 
of these fc-body gate can be approximated to accuracy e 
by a product of at most log*^ i one- and two-qubit gates 
taken from a finite discrete set containing, say, M gates, 
and where c is a constant (close to 2). 

Thus, to count the number of states that can be pro- 
duced by arbitrary time-dependent Hamiltonians, it suf- 
fices to count the number of polynomial-size quantum cir- 
cuits constructed from a universal discrete set of one- and 
two-qubit gates, and to consider an e-ball around the out- 
put of each of these circuits, i.e. the set of states within 
a distance e of the outputs of these circuits. Surely, 
the states reached in polynomial time by arbitrary time- 
dependent Hamiltonians are contained in the union of 
these balls. 

Since we limit the evolution to polynomial time, there 
exists a constant a such that the total number of gates in 
the simulation circuit is bounded by iiT" where K oc N 
is the number of qubits required for the simulation. 



There are no more than A/'circuits = (MK^)^ distinct 
ways of arranging these gates into a quantum circuit 
{M possibility for each gate and possible pairs of 
qubits between which it can be applied), and therefore 
no more than A'circuits distinct states that can be pro- 
duced. On the other hand, the states of ii'-qubits live 
on a (2^+^ — l)-dimensional hypersphere, whose surface 
area is 5 = 2tt^'^ /T{2-^), and an e-ball around a given 
state is a {2^^^ — 2)-dimensional hypersphere of volume 
V = 2n^'^-'^e^'^^'~^/r{2^). Combining, we see that the 
e-balls of physical states occupy only an exponentially 
small fraction ^""g'*"'"^ = 0{K^e^ ) of the total volume 
in Hilbert space. Thus, we conclude that the overwhelm- 
ing majority of states in the Hilbert space of a quantum 
many-body system can only be reached after a time scal- 
ing exponentially with the number of particles. For this 
reason, even for moderate size systems, these states are 
not physically accessible so they do not represent any 
state occurring in nature. For very large systems like in 
our expanding universe, it is indeed fair to assume that 
the time a system can evolve is not exponentially larger 
than the space it occupies. 

Conclusion — We have demonstrated that any time- 
dependent local Hamiltonian can efficiently be simulated 
using a quantum computer, independent of the frequen- 
cies involved. As an application, we have shown that the 
set of quantum states that can be reached from a product 
state with a polynomial time evolution of an arbitrary 
time-dependent quantum Hamiltonian is exponentially 
small. This means that the vast majority of quantum 
states in a many-body system are unphysical, as they 
cannot be reached in any reasonable time. As a conse- 
quence, all physical states live on a tiny submanifold, and 
that manifold can easily be parameterized by all poly- 
sized quantum circuits. This raises the question whether 
is makes sense to describe many-body quantum systems 
as vectors in a linear Hilbert space. The recent advances 
in real-space renormalization group methods [1, H^, Hfj 
indeed seem to suggest that a viable approach consists of 
parameterizing quantum many-body states using tensor 
networks and quantum circuits. 
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Appendix A: Average Hamiltonian 

Define the average Hamiltonian for time bin 

rAt 



where the last step folloows from Schrodinger's equation. 



H^^ ^— / H{s)ds 



(Al) 



and consider the time evolution operator C/^^(s) — 
exp{— isi?^^} that it generates. Our goal is to show that 
U^^{At) is close to [/(At) = rexp{-i H{s)ds} de- 
fined with a time-ordered integral. 

Let X{s) [/t(s)J7^^(s). We seek a bound on the 

norm of X{At). Since X{0) — 0, Schrodinger's equation 
yields 



X{At) 



At 



ds X{s) 



At 

I I dt [/t(s)AiJ(s)J7^^(s) 
'o 



with AH{s) = H{s) - H^^ . Using the definition 
Eq. (|Aip . this expression becomes 



\\U{s)-U{t)\\^ 



H{t)dt 



(A2) 



Appendix B: Decoupling principle 

Suppose that we replace each term Hx{t) in Eq. ^ 
by its local time average, that is, we define Hx{t) — 
J G{t — s)Hx{s)ds where G is a normal distribution of 
width a to be specified later, and we assume throughout 
that integrals are over the entire real axis unless other- 
wise stated. If T is the total evolution time, we extend 
the definition of H{t) so that H{t) = for t < or t > T. 
Note that these integrals involve only /c-body terms, so 
they can efficiently be carried-out classically Q. We sim- 
ilarly define H{t) = J2x ^^i^)- i*^®^ behind the 
convolution is to implement a low-pass filter with a fre- 
quency cutoff of order a~^. As we will now demonstrate, 
this has very little effect on the evolution of the system 
when a is sufficiently small. 

Formally, we will compare the time-evolution operators 
U and U generated by H{t) and H{t) respectively. Let 
X{t) = 1- W{-oo,t)U{-oo,t). We seek a bound on 
the norm of X{oo). Since X{~oo) — 0, Schrodinger's 
equation yields 

X{oo) = j dt X{t) = -i j dt U\t)AH{t)U{t) , (Bl) 

with AH{t) — H{t) — H{t). By commuting the order of 
the integrals in t and s in Eq. (IA2p . and using the fact 
that G integrates to 1 and is symmetric, we can bound 
||Ar(oo)|| from above by 



dt / dsG{t- s)\\U\t)H{t)l}{t)-U\s)H{t)lJ{s)\\ 



Using the triangle inequality for the operator norm, we 
obtain 



\X{oo)\\ < \\H\\ j dt j ds G{t - s)x 
x[\\U\t)-W{s)\\ + \\U{t)-U{s)\\], 



At I- At where = sup^ ||i?(t)|| and we used the fact that H{t) 

'IT I I {U\s)H{s)U^^{s) - U\s)H{t)U^^{s)) dsdiis supported in [0, T]. From Schrodinger's equation, it is 
Jo Jo easy to show that WUUt) ~ UUs)\\ < \t - s\\\H\\ and 



We can change the order of the integrals in the first term 
of this equation and use the triangle inequality of the 
operator norm, to conclude that 



easy to show that \\W(t) ~ U''{s)\\ 
U{t)-U{s)\\ < \t-s\\\H\\ < |t-s|||iJ||, where the last 



\XiAt)\\ < 



\H\\ 



At pAt 



At JO 
+ ||C7A^(.)- 
< 2\\HfAt 



{\\U{s)-U{t) 

U^^{t)\\) dsdT 



^ Although we want the Hamiltonian H(t) to be as general as 
possible, we need to assume that these integrals exist. Piece- 
wise continuity is enough to ensure this. If these integrals do not 
exist only due to some short time intervals of total duration e, 
then these intervals can be left out of the simulation at the cost 
of an error at most 
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step uses Young's inequality. Combining these bounds, 
we obtain 

||X(oo)|| <2\\Hf dt j ds G{t - s)\t - s\ 

<2\\HfT^a . (B2) 

This is the desired result: choosing the width a of the 
normal distribution to be much smaller than T||iJ|p im- 
plies that the norm of X{oo) — the difference between the 
exact evolution operator and the one obtained after trun- 
cating the high frequency components of H{t) — will be 
negligible. 



The bound in Eq. (jB2p is tight. To see this, con- 
sider the well known example of a two-level system with 
Hamiltonian H{t) = A[e*"*o-" + e"*'^*CT+], where cr" 
(ct^) is the lowering (raising) Pauli operator. When 
w ^ A and loT = 2Trk, for fc G Z, the evolution 
U{0,T) is well approximated by exp{iTX'^az/cj)- This 
is the so-called AC Stark shift. Averaging this Hamil- 
tonian over a window of width cr > — will truncate the 
unique oscillating terms, yielding H{t) « 0. Therefore, 
U{0,T) - C/(0,T) « expiiTX^a^/uj) - 1, and for short 
times T the total approximation error is of order X^T/lu, 
or ||_ff|pTcr, as claimed. This example and our general 
bound Eq. (|B2I) reflect the fact that high-frequency terms 
can only affect the overall evolution at second order in 
perturbation theory. 



